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The quark condensate is calculated in terms of the effective string tension and 
the constituent quark mass. For 3 colors and 2 light flavors, the constituent mass is 
bounded from below by the value of 460 MeV. This value is only accessible when the 
string tension decreases linearly with the Schwinger proper time. For this reason, the 
Hausdorff dimension of a light-quark trajectory is equal to 4, indicating that these 
trajectories are similar to branched polymers, which can describe a weak first-order 
deconfinement phase transition in SU(3) Yang-Mills theory. Using this indication, 
we develop a gluon-chain model based on such trajectories. 



I. INTRODUCTION 



At zero temperature, the classical action of QCD with A^f massless quark flavors possesses the 
so-called chiral SUL(Af)xSUR(Af) symmetry. On the quantum level, this symmetry is spon- 
taneously broken, and the order parameter for the symmetry breaking is the chiral conden- 
sate i^Tpip')- Together with confinement, characterized by the gluon condensate (^g'^{F^^Y^, 
chiral-symmetry breaking is one of the two most fundamental nonperturbative phenomena 
in QCD. We notice that the problem of interrelation between the chiral-symmetry breaking 
and confinement was already considered in Ref. [l]. In those papers, within a Nambu-Jona- 
Lasinio-type model with confinement incorporated, a proportionality of the two condensates 
was found. Here, we revisit this problem in a different framework. To this end, we use the 
definition of the quark condensate iipip') in terms of the one-loop effective action. 




By virtue of the world-line representation for r[Aj^] (cf. Ref. |2|), such an approach allows 
us to relate the condensation of a quark to the zigzag behavior of its trajectories. Quark 
condensation requires a I/a/s asymptotic behavior of the integral over the trajectories at 
large Schwinger proper times s [sl. To what extent this zigzagness is needed to provide quark 
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condensation and how to reconcile it with quark confinement are the questions addressed in 
this paper. 

n 

Being motivated by the above-mentioned proportionahty [1] between the chiral and the 
gluon condensates, we calculate (r[A^]) by imposing a confining area-law behavior of the 
quark's Wilson loop. To this end, we introduce a parametrization of the minimal area 
rendering the path integral, which represents ^r[y4j^]), analytically calculable. Then, the 
asymptotic behavior of this path integral in the chiral limit, together with the known 
expression for i^'ipip^ in the heavy-quark limit, leads to the 1/s-scaling of the effective string 
tension. We will show that this scaling corresponds to quark trajectories with the Hausdorff 
dimension equal to 4. Such trajectories, which are crumpled more severely than the ordinary 
Brownian random walk, are characteristic of the so-called branched polymers. To illustrate 
the physical significance of this class of trajectories, we have developed a model of the 
deconfinement phase transition, which is based on a gluon chain extending along such a 
trajectory. It turns out that this model describes a weak first-order phase transition, unlike 
its counterpart based on Brownian random walks, which describes a second-order phase 
transition [4|. As such, it is an interesting example of a QCD-inspired model, which can 
analytically describe the deconfinement phase transition in SU(3) Yang-Mills theory. 

The outline of the paper is as follows. In the next section, we develop a parametrization 
of the minimal area in terms of the contour of a Wilson loop and apply it to the calculation 
of the quark condensate. Starting with the purely exponential dependence of a Wilson 
loop on the minimal area, we then proceed, motivated by the heavy-quark limit, to more 
complicated ansatze. These ansatze, involving a pre-exponential area- dependence, allow 
us to recover within the same formalism the "area-squared" law of small Wilson loops. In 
this way, we find, within the present approach, the lowest possible bound on the value of 
the constituent quark mass, and demonstrate, both analytically and numerically, the 1/s- 
scaling of the effective string tension. In section III, we argue that this scaling corresponds 
to quark trajectories with Hausdorff dimension 4, and use such trajectories to develop a 
model of the deconfinement phase transition in SU(3) Yang-Mills theory. In section IV, we 
summarize the main results of the paper. In Appendix A, we check that the known result 
for the heavy-quark condensate is correctly reproduced by small Wilson loops, which depend 
quadratically on the minimal area. We furthermore show there that the ansatz for the Wilson 
loop involving the pre-exponential area- dependence can be used to fit the quadratic area- 
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dependence of small Wilson loops. In Appendix B, we illustrate for the simpler case of heavy 
quarks, how to get, within our approach, an equivalent confining Nambu-Jona-Lasinio-type 
model. 



II. QUARK CONDENSATE FROM THE MINIMAL- AREA LAW 



We use Eq. (JT]), where the averaged one- loop effective action reads 

ds _2 



X 



{2S + l)Ni 
tr V exp 



(2) 



drT^ (A^i, - 

In this formula, P and A stand, respectively, for the periodic and the antiperiodic boundary 
conditions, so that Jp = J^^^^^^^^^^y j\ = The trajectories z^{t) obey the 

equation 

/ drz^{r) = {). (3) 

Equation ([3]) means that the center of a trajectory is the origin, i.e. the factor of volume 
associated with the translation of a trajectory as a whole is already divided out. In Eq. ([2]), 

is a generator of the SU(A'"c)-group in the fundamental representation, 5* is the spin of the 
fermion, which propagates along the loop (i.e., for a quark, S = 1/2), and A^f is the number 
of light-quark flavors. Furthermore, since the quark condensation can occur only due to the 
gauge fields, we have subtracted the free part of the effective action, so that (r[0]) = 0. 

The field-strength tensor entering the spin term ~ ip^ip^Fj^^^ can be recovered by means 
of the area-derivative operator ^j— acting on the Wilson loop jsl- This fact allows us to 
reduce the gauge-field dependence of Eq. ([2]) to that of the Wilson loop, 



{W[ 



tiV exp [ig / drT^A^^z^ 



as follows: 







X < exp 



dTipf^ip^ 



Vz„ 



6 



SSf,^{z{T)) 



{W[z,])-NA. 



(4) 



We proceed now with the construction of the minimal area of the Wilson loop (Vrf-^^]) in a 
such a way as to render the path integral analytically calculable. 
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A. A simple ansatz for the minimal area of the Wilson loop 

We ignore perturbative interactions of the quarks, which in any case do not affect the 
nonperturbative quark condensate. Hence, let us start by introducing a simple ansatz for 
the Wilson loop, such as the one given by minimal-area law: 

{W[z^]) = N, ■ e-''^'>^^'-^'^^\ (5) 

Here, a{s) is some proper-time dependent string tension of a quark, and Smini-z^t] is the area 
of the minimal surface bounded by the quark trajectory z^{t). In the purely spatial 3- 
dimensional case, the minimal surface can naturally be built out of infinitely thin segments, 
so that its area reads S^d = \ dT\7, x z| [6]. Clearly, this is just an area of the surface 
formed by a rotating rod of a variable length. As such, this area is minimal from among the 
areas of all possible surfaces which can be formed by rotating curved lines. Furthermore, a 
4-dimensional generalization of is given by 



1 r 

>S'4d = —7= / dT\e^y\pZxZf,\ 
2^v 2, Jo 



The constant in front of the integral in this formula is fixed in such a way as to have 84^ 
going over to 83^ when Z4, — Z4. = 0. Notice also that, owing to condition ([3]), the surface 
pieces forming the macroscopic surface merge at the origin. For this reason, from among all 
the areas of cone-shaped surfaces with the boundary Zfj_{T), the minimal one is automatically 
given by 5'4d. 

The functional S'4d[-2^] can furthermore be minorated as follows. Denoting 
we have 



/ dra^u =. 2 • / dra^^] < ■ / dra^^ < V2 dT\ \ap„\ j . 

-^0 \ f.<uVO J Jo Jo \p^, J 

n j n 

Next, applying the Cauchy-Schwarz inequality in the form ^ |aj| < \ n-^a^^ where in 
our case n = 6, we get 

V2 f drlj^W^M <^ [ dr J2< = '^f drW,.\- 
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Therefore, 



(6) 



In what follows, we assume that S*^™ is provided by the right-hand side of this inequality. 
Denoting the integrated surface element and the rescaled string tension respectively as 



E,,,, = / dra,,,, and cr = ais) 



2^12' 
we can write the Wilson loop, Eq. ([5]), in the simple form 

(l^[^^]) = iV,-e-^l^-'. (7) 

We will show that the I/a/s asymptotic behavior at large s of the path integral in Eq. ([2]) 
in the chiral limit leads to a calculable s-dependence of the effective string tension a. 

By means of an auxiliary integration, we can further proceed from to in the 

exponent in Eq. ([7]): 

{W[z,])=N, 







We now apply the Hubbard-Stratonovich trick, which amounts to introducing an auxiliary 
constant Abelian field B^^y-. 

= f . f n rdB^y] e-*-t^---. (8) 



oo 



The trajectory-dependence appears now in the form e~5^'"'^'"', which can be recognized as 
a Wilson loop of an electron of charge (—1) in the electromagnetic field with the strength 
tensor 

*^ ui/ — ^ fii/XpB Xp- (9) 
The A-integration can be performed analytically. 



d\\'/\ 







1 _|_ iliiii 



to arrive at the following Abelian- like form for the Wilson loop: 

- y ■ n 7-^- (10) 
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B. Calculation of the quark condensate 



We are left to calculate the effective action and the quark condensate ([T]), by using the 
obtained parametrization of the Wilson loop, Eq. ( JTOj) . To this end, we notice that the term 
B^u^^u = J^fiu Jq drz^Zu, in the exponent of Eq. (ITOi) . can be written with the help of the 
(Abelian) Stokes' theorem as 



where the integration on the right-hand side goes over an arbitrary surface bounded by the 
contour z^{t). Therefore, when plugging Eq. (fTOj) into Eq. (j4j), we have 



exp 



-2 



6 



e 2 



2 ^(lU^ flL/ 



exp 



Accordingly, the path integral in Eq. (jl]) becomes that of an electron of charge (—1) in a 
constant electromagnetic field: 



1 



Vz,, 



p J A 



-2 / dTijjf.^p^ 



Ss^^{z{t))_ 



e 2 



^ . 1 • i 



{Airs)' 



(11) 



The result of integration is given by the so-called Euler-Heisenberg-Schwinger Lagrangian Q 
(for reviews see (j, ISl), so that 



Eq. (11) 



[s'^abcot{as) coth(6s) — l] . 



(12) 



In this formula, 



and 



- + v/(E2 - H2)2 + 4(E ■ H)2 



H2 



E^ + v/(E2 - H2)2 + 4(E • H)2 



are related to the two invariants of the auxiliary electromagnetic field J^^j, as a? — IP' = 
E2 — H2, 02^2 = (E ■ H)2. Equation (fT2|) is a compact-form representation of the fermionic 
determinant, which can be expanded in the number of external lines of the electromagnetic 
field. Retaining in such an expansion the leading term, which corresponds to the one-loop 
diagram with only two external lines, we have 



r o2 



Eq. (11) 



(47rs) 



-a') + 0{s'\J',,\') 



(47rs) 
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Returning to the i?^j,-field, which is related to the J-'^j.-field through Eq. ([9]), and neglecting 
the terms 0{s'^\B^X), we can write 



Eq- (11) -i^E^i- 



This expression, along with Eq. (fTOj) . yields for the one-loop effective action, Eq. 



1 15 1 



127r2 8 (27r) 



ds e-"" ' 



s a"" 



n / 



B 



1 + ^ 



7/2 ■ 



In this formula, i? is a 6-vector, 



B = {Bi2, Bi3, Bi4, B23, B24, so that Bf^^ = B"^. 

Furthermore, the neglection of the terms 0(^s^\B^^\'^^ leads to the corresponding restriction 
of the integration range. 

Using now the known value of the solid angle in 6 dimensions, = tt^, and fixing Nc = 3, 
we arrive at a remarkably simple expression: 



15Nr 1 ds e 



(An) 



m?s rl/s 737 

-— / dB- 



x7/2- 



Accordingly, the quark condensate reads 

15iVf 



647r2 



■ m ds 



00 rn?s rl/s r>7 



Integration over B can be performed analytically, and we obtain 

3N{ 



4772 



■m / dse-"" 'a' ■ f[A{s)] 



(13) 



where 



A{s) 



2aH 



1 and -(5^ + 4^ + 6) 



We arrive now at the crucial observation [sj] that, in the small-mass limit, (V'V') can only be 
finite if 



3/2 



0-^ ■ /[A] -> ^ at s ->cx). 



(14) 
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where ctq is some constant parameter to be determined. Because of the relation o"^ = 2X^7 
we conclude that, at large s, the following equality should hold: 

^ = x, where x = 2{aos)^/^. (15) 

_ 1 



Now changing s — )■ x in o" = 2^^, we obtain 



a 



21/3 

(16) 



ao V ^^^^ ■ ^(3^) ' 

where A{x) is a root of Eq. (fTSl) . The function (^[A] = is plotted in Fig.[TJ It vanishes at 
A = and A 00, and has a maximum at A ~ 15.36 = A: G[A] ^ 0.51 = Gmax- Therefore, 
for X < Gmax, equation G[A] = x has two roots, Ai{x) G {0,A) and ^2(2;) G (A, 00), while 
for X > Gmax it has no roots at all. Thus, we consider only 

and plot in Fig. [2] the ratio (t/cto, Eq- (H^ . as a function of x, for the physical root Ai{x). 
It turns out that the root A2{x) decreases with x so rapidly that the product x'*/'^yl2(x) 
decreases as well. As a result, a/do increases with x, i.e. with s, which is unacceptable from 
the physical viewpoint. Instead, the root Ai{x) increases with x, and so does the product 
x^/^yli(x), yielding a physically acceptable decrease of the (5"/cro)-ratio with s. This ratio, 
corresponding to Ai{x), can well be fitted by the following function: 



a 



0-49 ^-i. ^ ^ tjmax 

,.0.99 X ^ „ , _s 

2 ^yj^ 



.'^0/ fit I at ^iffi^ < X < Gmax- 

The result of the fit is plotted in Fig. [31 Therefore, at large s of interest, we obtain the 
following scaling behavior: 

In the next subsection, we will show that, for a Wilson-loop ansatz more realistic than that 
of Eq. ([7]), the fall-off of a with s becomes linear and fit-independent. 

The existence of the maximal value Gmax of leads to a lower bound for the value of 
the constituent quark mass m. Indeed, from Eqs. (fT3|) and (fl^ . we obtain for the chiral 
condensate: 

(^^> = -^-ao I ^e --^-0 for -W = -^^J >1- 



0.6 




A 

Figure 1: The function G[A] = ^ for A G (0, 100). 




0.55 



X 

Figure 2: The ratio crja^^ Eq. (I16p . for the root A\{x) of the equation G\A\ = x. 
Equating further this expression to the commonly accepted phenomenological value, 

(^/'V) = -A/", M = (250 MeV)^ 
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The same as in Fig. 2 
The fitting function 
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Figure 3: The ratio a/aQ, Fig. [21 and the fitting function, Eq. pT|) . 



we get 



(To = TT 



2/3 



and m > 



Af 



1/3 



(19) 



3iVfG 

max , 

In particular, for = 2 light flavors, the lower bound for the constituent quark mass is 
m > 610 MeV. That is too big. The large value of this bound indicates that the purely 
exponential fall-off of the Wilson loop with the minimal area, assumed in Eq. ([7]), is too 
naive and should be improved. Moreover, the purely exponential fall-off of the Wilson loop 
with the minimal area is incompatible with the Gaussian fall-off arising for Wilson loops 
with sizes smaller than the vacuum correlation length jol (cf. Appendix A). It turns out 
that these two problems can be simultaneously improved. That will be done in the next 
subsection. 



C. Diminishing the lower bound for the constituent quark mass 

A natural way to improve on the purely exponential fall-off of the Wilson loop with the mini- 
mal area is to consider its additional pre-exponential area-dependence. Such a pre-exponent 
would not change the static quark- ant iquark potential. As will be shown in Appendix A, a 
suitable parametrization of a Wilson loop, which, for a constant string tension, accommo- 
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dates the area dependence for both small and large loops is given by the formula 

(^M> = 2^^^V{a) ■ (^l^^'^l)" ■ ^"(^I^M^I)' (20) 

where T{x) and Ka{x) stand, respectively, for Gamma- and MacDonald functions, and 
a > is some parameter. Equation f l20|) obeys the normalization condition (W^[0]) = N^., 
and reproduces the purely exponential fall-off, Eq. (I7l), for a = \. In Appendix A, we show 
that, for a constant string tension, the best fit of the quadratic area-dependence is obtained 
for a in the narrow interval a G (1.83, 1.9). In this subsection, by recalculating the chiral 
condensate with the s-dependent string tension instead of a constant string tension, we 
show that these values of a lead also to the lowering of the constituent quark mass down to 
460 MeV. 

The calculation of the condensate can be done by using for Eq. fl2Up an integral represen- 
tation 

AT 2 

and introducing the field B^^, through the Hubbard-Stratonovich transformation ([8]). We 
can then again analytically integrate over A, and obtain [cf. Eq. ffTOj) ]: 

We proceed along the lines of the previous subsection to obtain the following generalization 
of Eq. (DSI): 

{^^) = -m f^ds e-^'^ ■ ]\A{s), a], 

where 







~. , _ 4{6A(1 + A)"(2 + A) + 6[(1 + AY - 1] - (2 + a)A[6 + (1 + a)A(3 + aA)\} 
'"^ ~ 3(a-l)(l + 

In particular, one can check that /[A, 1/2] = /[A], and that, despite the [a — l)-term in the 
denominator, / has a finite limit a — 1: 

, 4[6(l + A)3ln(l + A)-A(llA2 + 15A + 6)] ^ 
j\A,a\^^ 3(/^ at «^1, 

i.e. / is a continuous function of a. In Fig. HI we plot the function (^[A, a] = /[A,^]/^ 
for three different values of the parameter a: a = 0.5 (in which case it coincides with the 
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function G[A] from Fig. [T]), a = 1, and a = 1.5. We observe a sharpening of the peak with 
increasing a. The highest peak can be found numerically to be: 

Gina.x[A,a] 1.18 at a > 1. 

The corresponding lower bound for the constituent quark mass, defined by Eq. (fT9l) with 
Gmax replaced by G^ax, is 

m^in = 460 MeV. (22) 

Furthermore, the equation G[A, a] = x again has two roots, one of which, Ai{x, a), leads to 
the decrease of a with the increase of x. Because of the sharpening of the peak of (^[A, a] 
with increasing a, this root behaves as Ai(x,a) ~ x'^^°'\ where e{a) — at a — > oo. The 
corresponding x-dependence of the effective string tension a is determined by Eq. (fT6l) with 
A(x) replaced by Ai(x, a — t- oo), and reads a ~ x~'^^^. Recalling the definition of x in terms 
of s, Eq. ( ITSl) . we conclude that 

a =- at a > 1, (23) 
s 

where c is constant. Thus, the proper-time scaling of a becomes strictly 1/s for a ^ 1. 
In practice, as one can see from Fig. HI this scaling holds already for a > 1, in particular 
for the values a G (1.83,1.9), which for the constant string tension best fit the quadratic 
area-dependence of small Wilson loops (cf. Appendix A). 

To find the value of the coefficient c in Eq. (p3l) . we notice that in the real world the string 
tension (jf (where "f" stands for fundamental representation) is equal to (jf = (440 MeV)^ 
at distances R < Rg.h., and vanishes at distances R > Rs.h.- Here, Rs.h. = ^ (with "s.b." 
standing for string-breaking) is the distance at which the string joining the quarks breaks 
by converting its energy to the production of another quark-antiquark pair, of mass 2m. 
Our model replaces such a step-function behavior of the string tension by the function a. 
Therefore, the normalization coefficient c can be found by considering the work necessary 
to adiabatically separate, during the proper time s, a quark from an antiquark, from the 
minimal distance a needed for the string to exist to the maximal distance i?s.b. which can be 
reached without string-breaking. While in the real world this work is equal to af{Rs.h.—ci), in 
our model it is J^""'^' dRa{R), and c can be found from the equality of these two expressions. 
To calculate the last integral, we notice that, during the proper time s, a quark going away 
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Figure 4: The function G[A,a] = for A G (0,100) and a = 0.5 (cf. Fig. ^, a = 1, and 

a = 1.5. 

from an antiquark separates from it by the distance R = ms. Therefore, cf[R) = c^, and 
we obtain 

2m — CTftt 
m ■ In (^) 

" n 

For the value a ~ 0.34 fm in SU(3)-theory with quarks [10|, and m = 460 MeV, Eq. (122|) . we 

get an estimate 

c ~ 1.26. 

Finally, we compare with Rg.h. the maximal distance msmax achievable by a condensed 
quark in our model. The value of Smax can be obtained from the formula x = 2(o"os)'^/^, 
Eq. (fT5|) . where ao is given by Eq. (fT9l) . and the maximal value of the variable x is either 
Cmax or Gjnax; depending on the ansatz for the Wilson loop. For instance, for Xmax = Cmax, 
we have 

1 fSNiG 

max 



'""^ 47r V AT 

Furthermore, because of the common factor of m in both Rs.h. = 2m/ ai and mSmax, it 
amounts to comparing 2/crf and Smax- We obtain 

2/o"f I 3.85 for Xmax = Cmax, 
■Smax I 2.20 for a;max = Gmax- 
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Thus, at the distances which a quark in our model can achieve, the string does not break. 
We proceed now to a discussion of the physical meaning and applications of the obtained 
scaling of the effective string tension, Eq. (l23|l . 



III. HAUSDORFF DIMENSION OF A LIGHT-QUARK TRAJECTORY. A 
MODEL OF THE DECONFINEMENT PHASE TRANSITION IN SU(3) 

YANG MILLS THEORY. 

In this section, we argue that the linear fall-off with large proper times s of the effective 
string tension is associated with a highly fractal behavior of quark trajectories at large 
distances, and make quantitative predictions out of this behavior. In a constant electric 
field E, classical Euclidean trajectories of a particle of mass m and electric charge (—1) are 
circles of radius -Rjnst = rn/E. (The subscript "inst" is due to the fact that these classical 
trajectories are commonly referred to as world- line instantons ll|.) The path integral for 
the leading fluctuating trajectory about the world-line instanton is Gaussian [11]. Therefore, 
the Hausdorff dimension of such a trajectory is the same as for a free Brownian random walk, 
i.e., 2. That is, L ~ -Rfnsf where L is the length of a fluctuating trajectory. In our case, 
E a 1/s, therefore -Rinst ~ f^s, as a direct consequence of the linear fall-off of the 
effective string tension a with s. Thus, we have 

L ~ (ms)^ (24) 

Let us now take a quark trajectory passing through two different points, x and y, and 
consider a free random walk from x to y, which takes place along the trajectory during the 
proper time s. The length L^y of this random walk is of the order of L, while the end-to-end 
distance Rxy = \x — y\ for a free random walk scales with the proper time as R^y ~ s^^^. 
Therefore, we have from Eq. (12^ : 

^xy ~ Ll^'y^, 

i.e. the Hausdorff dimension of the fluctuating trajectory of a light quark is 4 instead of 2. 
We emphasize once again that this result is a consequence of the linear fall-off with proper 
time s of the effective string tension 5" of a light quark. 

The behavior L^y ~ R'^y indicates that trajectories of lighter quarks are much more 
crumpled than the conventional Brownian random walks, for which L^y ~ R^y. Specifically, 
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trajectories with Hausdorff dimension 4 are characteristic of the so-called branched poly- 
mers 12|. These are just polymers whose monomer branches are Brownian random walks 
themselves l^- Let us study the case of pure SU(3) Yang-Mills theory, where the decon- 
finement phase transition is weak first-order. We will show that such a phase transition can 
be described by a gluon chain connecting static quarks [l^, when this chain has Hausdorff 
dimension 4. In particular, we will reproduce the value of the critical exponent 

V = 1/3, 

which characterizes a weak first-order phase transition |l5|. 

To this end, we follow the approach of Ref. j4|. In that paper, a gluon-chain model of 
the quark- ant iquark string, based on the Brownian random walk, was proven useful for a 
description of the second-order deconfinement phase transition. Here, we generalize this 
model for the case of the random walk with Hausdorff dimension 4. It amounts to noticing 
that the length of a random walk is L = a- (number of steps), where in the present case 
(number of steps) = ^. (In these and subsequent formulae, a is the length of one step, i.e. a 
mean distance between the nearest gluons in the chain.) Thus, 

Next, since each segment of the chain between the nearest gluons may have its own color, 
the entropy of the whole chain is = \n(^Nc^°'^, where L/a is the total number of segments. 
The free energy of the chain reads -Fchain = <^iL — TS. The partition function of the random 
walk has the form 



PC 

„ Jo 



ds 



exp 



+ (/3n) 
4s 



chain 



,0 (4vrs) 

where the first term in the exponent corresponds to the kinetic energy of the walker. Here, 
i? is a 3-dimensional quark-antiquark separation, /3 = 1/T, and the sum over winding modes 
n runs from — oo to -|-oo. This partition function can therefore be written as 



ds 



(Airs) 



exp 



+ {/3n) 
4s 



InNr 



Changing the variable s — )■ 1/s and carrying out the resulting saddle-point integration, we 
obtain 



1 



Pert 



\nNr 



1/3 



X 
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X 



[R2 + (/3n)2]2/3 



■ exp 



1/3 



Ir2 + (/3n)Y' 



Aa \ a 

At large quark- antiquark separations R of interest, only the zeroth winding mode survives, 
so that 



1 



1 



4 V 37r3a 



PCTf 



i?4/3 



exp 



4a 



InNr 



1/3 



(25) 



at i? — 7- oo. As follows from the exponent in this formula, our random- walk model describes 

I ! 

a static potential, which grows at large distances as R^/^ [20]. That is, Vt{R) = 'y{T)R^/^, 
where the dimensionality of the coefficient 7(T) is (mass)''/'^. The critical behavior of 7(T) 
defines the sought value of u. Following the approach of Ref. [4(|, we deduce this critical be- 
havior from the full free energy of the system. The latter is the sum of the zero-temperature 
potential, 

Vt=o{R) = lR^'\ (26) 
and the free energy of the random walk, so that 

/rr\ ^ 1 Zr.^{R,T) 

7(T) = 7 - ^ttttIu 



i?— >-oo 



i?4/3 Z,„.(i?,To) 

The free energy of the random walk is normalized by the condition Z^.^ {R, Tq) = 1, where 
To is the temperature below which the random walk is suppressed due to the weakness of 
string fluctuations. (In what follows, we will prove that Tq is by a significant factor of 4.4 
smaller than the deconfinement critical temperature T^.) From Eq. (l25l) we have 



7(T) = 7 + 



3T 
4^ 



1/3 



- /3oo-f 



1/3 



(27) 



thus obtaining the deconfinement critical temperature, which happens to be the same as for 
the second-order phase transition 

a{a 



The temperature Tq is fixed by the condition 7(Tc) = 0, yielding 

To 



1 



(Tf 



IniV, 



^ ' 3(jf 



(28) 



(29) 



(InA^e 



Thus, we find the following critical behavior: 



7(T) ^ - 



at T — Tc, 
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i.e. the critical exponent v is equal to 1/3, that corresponds to a weak first-order phase 



transition 



15|- 



Furthermore, as follows from Eqs. (!28|) and ( l29l) . the temperature exceeds the temper- 
ature To by a factor of 

^ = l + a-f^^V(lniV.r. (30) 



To get a numerical value of this expression, we evaluate 7 in terms of o"f. This can be done by 
approximating the zero-temperature potential (l26l) as Vt=o(-R) — 7 (-R^^^) ■ -R, and equating 
it to the standard linear potential oiR. Here, we have used the obvious fact that B)-!^ varies 
much slower than R. The mean value (^R^^^'j corresponding to the potential f l26p reads 

^ irdRR'/'e.p[-^R-/^.T] 
^ ' j;^dRexp[-jR^/^-r] ' 

where T is the Euclidean time needed to communicate a signal between the two static 
quarks. For our evaluation purposes, it suffices to set this time to its minimal value, that is 
T = R. The corresponding mean value of R^^^ then reads 

/^i/3) ^ 3 ■ r(4/7) _ 

\ / " 7 -1(10/7) ^ ' 

where T(x) again stands for the Gamma-function. The matching condition, 7 (^R^^^'j = ai, 
yields the desired expression of 7 in terms of af. 

"7 -1(10/7) T'/' 



7 



CTf 



1.4 - a'/'. 



3 - r(4/7) 

Plugging it into Eq. ( !30|) and setting in that equation = 3, a = 1 GeV~^ lO, ll6| in SU(3) 
Yang-Mills theory, and o"f = (440 MeV)^, we obtain 



^ ^ 4.4. (31) 

-1-0 

Thus, the relative range of temperatures, in which our random-walk model of the deconfine- 
ment phase transition is applicable, is rather broad. 

Finally, let us consider the limiting case Nc = 1, i.e. an Abelian theory with confinement 
based on the condensation of magnetic monopoles. In such theories the role of gluons which 
the chain passes through is played by the off-diagonal gauge bosons. Upon the spontaneous 
breaking by an adjoint Higgs field of the initial SU(2) symmetry down to (compact) U(l), 



these bosons become massive 



17| . Since the original non- Abelian symmetry is SU(2), one 
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expects that the phase transition is second order. For a gluon-chain model based on the ran- 
dom walk with Hausdorff dimension 2, this is anyhow the case [4], since the phase transition 
in that model is second order for Nc > 1 as well. Instead, for the model developed above, 
where for Nc > 1 the phase transition is weak first order, the reduction to the second-order 
phase transition for A^^c = 1 is not a priori obvious. To demonstrate it, we write Eq. (1271) for 
Nc = l as 

=7 + ^^""^2^ ■x^(l-x), where x = (T/To)^/^ 
The critical value Xc is defined by the equation '~f{xc) = 0, which explicitly reads 

To derive the critical behavior of 7(x), we Taylor expand this function near Xc, by sub- 
stituting X = Xc — C,, where ^ ^ Xc- To the first order in ^, this expansion yields 
7(x) — ^ • Xc{3xc — 2). We note further that Xc{3xc — 2) = f'{xc), where the function 
f[x) = x'^{x — 1) is positive-definite together with its derivative for any a; > 1. Since 
Eq. (132|) . defining Xc, can be written as f{xc) = ^(^J^^l^i/s i which is also a positive-definite 
quantity, we conclude that Xc > 1, and hence f'{xc) > as well. Therefore, Xc{3xc — 2) > 0, 
and 

7(x) oc {xc — x) (X 

where the proportionality coefficients are positive-definite. Representing now T in the vicin- 
ity of Tc as T ^ Tc(l — y) , where r ^ Tc, we get 

3Tc' 

where the proportionality coefficient in front of (Tc — T) is positive-definite too. Therefore, 
we conclude that 

7(T)oc(T,-T), 

i.e. the phase transition is second order, with the critical exponent u equal to 1. We have 
thus demonstrated that the weak first-order phase transition, described by our model at 
Nc > 1, becomes for A^^^ = 1 second order, with the universality class of 2-dimensional 
Ising model. The second-order phase transition at A^c = 1 is due to the insufficient entropy 
generated by the random walk in that case, even for a walk with Hausdorff dimension 4. 
As a result of this insufficiency, the model from the thermodynamical viewpoint becomes 
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equivalent to the conventional model of the deconfinement phase transition based on the 
condensation of long cl.se,. tn.^ ^ The latte. also belongs to the 2d-Ising nn.ve.saUty 
class (for a discussion of this issue see |^). 



IV. SUMMARY 



It is intuitively clear that the condensation of quarks is associated with the zigzag character 
of their motion, leading to a decrease with the Schwinger proper time of their effective string 
tension. It is also known that small Wilson loops of heavy quarks respect the "area-squared" 
law, which reproduces the heavy-quark condensate. In this paper, using an ansatz for the 
Wilson loop interpolating between the area law at large distances and the "area-squared" law 
at small distances, we have shown that the decrease of the effective string tension with the 
Schwinger proper time occurs linearly. This leads to the lower bound ( l22l) for the constituent 
quark mass in our model, whereas the naive use at all distances of the area law alone yields 
an unacceptably high value of this mass. 

The linear fall-off of the effective string tension with the proper time is associated with 
highly fractal trajectories of light quarks. Quantitatively, we have found that these trajec- 
tories have Hausdorff dimension 4. That is, when such a trajectory connects two different 
points, the length of a random walk performed along this trajectory grows as the fourth 
power of the end-to-end distance. Trajectories of this kind are characteristic for the so- 
called branched polymers, each segment of whose is a Brownian random walk. We have 
furthermore used these specific trajectories to develop a gluon-chain model of the weak first- 
order deconfinement phase transition characterized by the critical exponent u = 1/3. Owing 
to the Z3-symmetry universality argument, supported by numerous lattice simulations, the 
deconfinement phase transition in SU(3) Yang-Mills theory is commonly accepted to be of 
this kind. For the realistic values of the string tension at and the vacuum correlation length 
a, we have demonstrated a relative broadness of the temperature range where our model is 
applicable [cf. Eq. (13T]) ]. Finally, in the limiting {Nc = l)-case, the deconfinement phase 
transition described by our random-walk model becomes second order, as it should due to the 
center-group universality argument [i.e. because the compact U(l) group emerges by means 
of a spontaneous breaking of an SU(2) symmetry, and this symmetry forces a second-order 
phase transition]. 
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In conclusion, our paper suggests a parametrization of the Wilson loop compatible with 
both the heavy-quark condensate at small distances and the chiral condensate at large 
distances. It provides a lower bound for the constituent quark mass, and quantifies fractal 
properties of the light-quark trajectories. Furthermore, the paper develops a model of the 
deconfinement phase transition in SU(3) Yang-Mills theory. It would be very interesting to 
explore how the obtained 1/s-scaling of the effective string tension changes with temperature, 
as this change should eventually lead to the chiral-symmetry restoration. 
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Appendix A. Fixing parameter a by the heavy-quark limit. 

In this Appendix, we illustrate how to fix the parameter a entering ansatz ( l20l) . To this end, 
we notice that this ansatz can fit the "area-squared" law, which describes the nonperturbative 
contribution to small Wilson loops jol. Let us first rederive this law and show how it 
reproduces the known expression (l9| for the heavy-quark condensate. 

By using the non-Abelian Stokes' theorem and the cumulant expansion, one can write 



the Wilson loop in the fundamental representation as ^| 



(iy[z^]) ~ tr exp 



(Al) 



where the integrations go over the minimal surface bounded by the quark trajectory 2;^(t), 
and ^xx' = V exp^ig J^^ T"y4^(M)(iu^] is the phase factor along the straight line. Further- 
more, the factor of 1/2! in Eq. (A.l) is due to the cumulant expansion, while the factor of 
1/4 is due to the (non-Abelian) Stokes' theorem. Now, for surfaces fitting into the circle of 
the size of the vacuum correlation length a, the field-strength tensor F'^^{x) can be treated 
as constant, and one approximates the correlation function in Eq. (A.l) as 

(F;,(a;)T'^$..,Fi(x')T^$.'.> ^ ■ ^ (5,a5.p - 5,Aa) ■ tr(W^) ■ {f;,{x)F%{x)) . 



21 



Using this parametrization in Eq. (A.l) and noticing that tr (T'^T ) = hS , one has 



{W[z^])-N,-exp 



96iV, 



{A.2) 



where (C^) = {g'^{F^yY) is the gluon condensate. In particular, for a flat non-selfintersecting 
trajectory ^/^(t), the double surface integral reads j do'^i,{x) J d(7^,^{x') = 23"^, where 5* = 
S[z Jis the area of the flat surface bounded by z^{t). Thus, one gets the "area-squared" 

law y 



{W[ 



48iV, 



{A.3) 



~ Nc ■ exp 

flat, non— selfintersecting Zij,(t) 

We apply now Eq. (A. 2) to the calculation of the heavy-quark condensate. To this end, 
we rewrite the double surface integral as 



d(r^y{x) / da^y{x') 



da^^{x) I da^p{x')d^d^p'{x - x') 



and use the (Abelian) Stokes' theorem to obtain 



- (b dz„ (b dz'iz - z' Y . 



doy,y{x) J daf,u{x') = 

One can notice that only the (z2;')-term in {z — z')"^ yields a nonvanishing contribution to 
the integral, so that 

2 



da^y{x) I da„Jx') = -- (p dz,y i dz[,(-2zz' 



dzi,Zp 



Inserting this expression into Eq. (A. 2), we can apply the Hubbard-Stratonovich trick, 
Eq. (IHl), to further write the Wilson loop as 



{WW]) 



n 



dB 



e Tec f ^MJ^ f dzuz^ 



(AA) 



where we have denoted for brevity C = -^gj^r- The Wilson loop in this form can now be 
used in the one-loop effective action, Eq. (|4]), to obtain 



Vz^ / Dip exp 



2NM 



d^ ^-m'^s 



n 

\fl<U 



+ 00 



dB 



e X 



1 1 ■ i 

d'T { ^zl + -'ipf.ipf, + -B^^z^Zu - iBf,^i:^i)i, 



(Aits) 



{A.5) 
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The last path integral coincides with that of Eq. (fTT]) . with the substitution J^^y — t- B^,^. 
Therefore, it is given by the Euler-Heisenberg-Schwinger Lagrangian, which to the order 
yields for the curly bracket in the last equation: 

-,2 



{■■■1^ ^ .-Yb^ 



3 



Applying now definition ([T]) and recalling that the solid angle in 6 dimensions is fig = vr 
we obtain the heavy-quark condensate in the form 

AN N( 1 TT^ /"°° 9 f^/'' ^ b2 



h^'^-y (87rC)3 (47r)2 3 

Because of the factor e""*^'', the essential values of the proper time are s < therefore 
7 > m^. In particular, in the leading large-m approximation, one can approximate the 
upper limit of l/s in the last integral by +oo, thus decoupling the s- and the 5-integrals 
from each other. The i?-integral then reads 



„2 

dBB'^e'^ = 12288 C^, 



leading to 



Finally, recalling that C = ^^j^' , and writing (C^) as (C^) = 47ra.s((F^",)^), where a,, = tz 
we recover the known expression for the heavy-quark condensate 



47r' 
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heavy 



-Nf 



a 



s \ 



127rm 

Afc=3 



This calculation proves that the "area-squared" law, Eq. (A. 3), for the nonperturbative 
contribution to a small Wilson loop is consistent with the large-mass limit of the quark 
condensate. 

To evaluate parameter a in Eq. f l20l) . we replace in that equation 5"|S^i,| — )■ a^S. Intro- 
ducing in the "area-squared" law, Eq. (A. 3), instead of 5 a new variable x = afS, we now 
want to find a, which provides the best approximation in the formula 




23 



for X < 1. We notice first of all that, up to the positive-definite coefficient 2"-iV(q;) ' 
second derivative of the approximating function reads 

— (x'^KJx)) = -2"-2r(a - 1) + a(2a - 1) • 2-"r(-a)s2(°-i) + Cfx^") + 0(x), a ^ 1. 

{A.8) 

Therefore, for a > 1, the second term on the right-hand side of Eq. (A.8) is subleading 
compared to the first term, which is negative-definite. Rather, for a < 1, the second term 
on the right-hand side of Eq. (A.8) becomes the dominant one, and one can check that the 
expression which defines its sign, that is (2a — l)r(— a), is negative-definite for a G (|, l), 
while becoming positive-definite for a < |. In the special case of a = 1, ^ {x°'Ka{x))\^_^ = 
Inx + 0{1), that is also negative-definite for x < 1. Therefore, we conclude that for any 
a > |, as considered in subsection II C, one has £^ (x"ir„(x)) < 0, in the same way as for 
the left-hand side of Eq. (A. 7). This finding suggests the simplest, analytic, way of seeking a 
by comparing the leading (9(x^)-terms on the two sides of Eq. (A. 7). Indeed, one can check 
that the (9(x^)-term on the right-hand side of Eq. (A. 7) is the dominant one for a > 1, in 
which case 

1 ^2 



-x'^K^ix) = 1 + + 0{x"^) + 0{x' 



2°-ir(a)"" " ■ 4(1 -a) 

Comparing this expansion with 



i^x^Vl ^x^ 



we find 



a = 1 -I- 



' 

that is indeed larger than 1. Substituting into this formula Nc = 3, (C^) = j^, 
where the value of the vacuum correlation length in QCD with dynamical quarks is a = 
1.72 GeV~^ (lOj, we obtain 

a ~ 1.90. 

A very close value 

a ~ 1.83 

can be found by numerically minimizing fo^' the difference between the two sides of 
Eq. (A. 7) in the relevant range x G (0,0.45). Here, the maximal value x^^x — 0.45 corre- 
sponds to the maximal radius a/2 of a circular contour for which the "area-squared" law is 
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Figure 5: The left-hand side of Eq. (A. 7) and its right-hand side for a = 1.9 and a = 1.83. 
still applicable. Indeed, the area encircled by such a contour is S = 7r(a/2)^, that yields 

a^max = CTlS ~ 0.45. 

In Fig. m we plot the left-hand side of Eq. (A. 7) and its right-hand side for a = 1.9 and 
a = 1.83, and observe a good agreement between these three curves. We conclude by noting 
that, for a G (1.83, 1.9), the ansatz 

approximates the "area-squared" law rather well. At the same time, as we have seen from 
Fig. m the linear fall-off of a with s holds with a high accuracy for any a > 1, in particular 
for a's from the above-mentioned range. 



Appendix B. A relation to the effective quark models. 

In this Appendix we derive, as an illustration, the Nambu-Jona-Lasinio action for the simpler 
case of heavy quarks. To this end, we represent the exponentiated Eq. (A. 5) before the -8^,^- 
averaging as 

expl-2N^NiJ —e-'^'x 
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X 



J A 



exp 



- / dr 



1-2.1^ : . ^ 



det(7^£)^ + m) 
det(7^a^) 



(5.1) 



Here, the Grassmann fields ip"" and ip'^ represent auxiliary Abelian spin-1/2 fermions. In order 
to represent the number of quark degrees of freedom, these fields are supplied with additional 
quantum numbers a — 1, ... , N^Nf. The covariant derivative has the form = 9^ + iC*^, 
where C^{x) — ^x,^B^n is the vector-potential corresponding to the x-independent field- 
strength tensor 3^,^. We apply now to both sides of Eq. (B.l) the 5^,^-averaging defined 
as 



n 



+ 00 



e 



Since B'ij^ = 2 ^ S^^, we see that this average has the correct normalization, i.e. (1) 



Furthermore, owing to Eq. (A. 5), the one- loop approximation adopted throughout the paper 
yields, for the left-hand side of Eq. (B.l), the expected expression exp ((r[A^])) that is the 
gluonic-field dependent part of the one-loop QCD partition function. On the right-hand side 
of Eq. (B.l), we have 



n 

\H<v 



+00 



dB 



e 16C X 



X 



The second exponent in this expression, which describes the subtracted free part, is of no 
relevance to chiral-symmetry breaking. The S^i^-average of the first exponent yields the 
desired action of a Nambu-Jona-Lasinio-type effective quark model. We note that the 
Sjuj,- integration is the same as in Eq. (A.4), up to the substitution 



^ dz^Z/^ J 



d^x x^'il)°--i^'il)°-. 



{B.2) 



For this reason, we can benefit from the already known fact that Eq. (A.4) is just a rep- 
resentation of Eq. (A. 2). This leads to the following action of a Nambu-Jona-Lasinio-type 
model: 



/ 



96Nr 



d'x / d^'x' x^x'^irxi.rxm'i.rx')- 
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The direct product 7j, ® can further be disentangled by using the standard Fierz 
transformation. What is important for the present discussion is the quadratic coordinate 
dependence of the four-fermion interaction. It has been shown in the first paper of Ref. jl| 
that, in the chiral hmit, this type of interaction leads to the chiral condensate i^ijjil)^ oc 
—a ■ ((7^E^), where E is the chromo-electric field. The reason why it is the chromo-electric 
rather than the full condensate entering the result is the so-called modified coordinate gauge 
adopted in in which ^4(0, X4) = Xiy4j(x, x^) = 0. Fixing this gauge, one can eliminate the 
chromo-magnetic part of the full condensate. That part is unlikely to be related to chiral- 
symmetry breaking, since at the deconfinement critical temperature the chiral symmetry is 
restored, whereas the chromo-magnetic condensate not only survives but can even increase 
further with temperature. At zero temperature, relaxing such a gauge-fixing, one obtains 
the chiral condensate (^ipip') oc —a ■ (C^). In this formula, the vacuum correlation length a 
plays the role of the inverse UV cut-off, which in the heavy-quark limit at issue is replaced 
by 1/m. That eventually leads to the QCD sum-rules result for the heavy-quark condensate, 
Eq. (A. 6). Remarkably, the A^f-proportionality of (^'^)heavy recovered by the sum over a 
in the Abelian condensate (^ip""ip"'^ . 
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